272 


0G5.4-1 


A Simulation of High Energy Cosmic Ray Propagation I 

* * ** * 
M. Honda, T.Kifune, Y. Matsubara , M. Mori , K. Nishijima and M. Teshima 


Institute for Cosmic Ray Research, University of Tokyo, Tokyo, 188 Japan 
*Department of Physics, Kyoto University, Kyoto, 606 Japan 

**The Graduate School of Science and Technology, Kobe University, Kobe, 657 
Japan 


Abstract 

^ p_We ^mulate high energy cosmic ray propagation of the energy region 
10 * -10 eV in the inter steller circumstances. In conclusion, the 

diffusion process by turbulent magnetic fields is classified into several 
region by ratio of the gyro-rgd^us and the scale of turbulence. When the 
ratio becomes larger then 10 ’ , the analysis with £he assumption of point 

scattering can be applied^w^th the mean free path E . However, when tfye 
ratio is smaller than 10 , we need a more complicated analysis or 

simulation. Assuming the turbulence scale of magnetic fields of the Galaxy 
is 10-30pc and the mean magnetic fields strength is 3 y^cgo gauss, the 
energy of cosmic ray with that gyro-radius is about 10 eV. 


1 Introduction 

Ij/I^y authors have studied cosmic ray propagation in magnetic 
fields' 1 . Especially, Gleeson and Axford had solved a diffusion equation 
with the assumption that cosmic ray are scattered, bv scattering centers. 

The recent observation of Galactic magnetic field ^ shows the average 

field strength is about 3 micro gauss and the turbulent fields are about 
1.5 micro gauss. The gyro-radius of cosmic ray in the energy region, in 
which we are interested, becomes 0.1-300pc. The random magnetic fields 
cannot deflect the cosmic rays so hardly, but an amount of small deflection 
may cause a large scattering angle as a result. We know that, in some 
cases, we can sum up the effects in terms of point scattering and mean- 
free-path(mfp) . However, we are not sure this is the case in the cosmic ray 
propagation in the turbulent magnetic fields, especially when the gyro- 
radius of the cosmic ray becomes near to the scale of irregularities. The 
assumption of point scattering can be an over simplification in the study 
of cosmic ray propagation in that energy region. It should be checked by 
the numerical calculation or simulation. 

In this paper, we present a Monte Carlo simulation of cosmic ray 
propagation in the Galactic space and show a result in terms of diffusion 
tensor. Gleeson and Axford showed the diffusion tensor have three 
independent element (K,,,Kj_, K T ) . The analytic form of these element as a 
function of mean free path are given by 
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where \ is the mfp and CO is the gyro frequency of cosmic rays in that 
magnetic field. When we leave from the idea of point scattering of cosmic 
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ray, we should not use mfp at the beginning of the study. Rather it should 
be given as the result of diffusion process. In our simulation, we 
calculate the two diagonal elements K . . and Kj^ . We check the consistency 
comparing these formula. When the consistency is maintained, we try to 
calculate the mfp from the formula. 


2 Generation of random magnetic fields 

It seems to be a difficulty in the simulation with irregular magnetic 
fields to make the flux line conserved. This difficulty is avoided in our 
case by generating the vector potentials by random numbers . The steps to 
generate the random magnetic fields are following. First we embed a lattice 
in the simulation space, whose lattice constant is chosen as the scale of 
the irregularity of the magnetic field. We generate the vector potential 
A=(A , A , A ) for each lattice point with the probability given by 

X Y Zj 

exp(-AVt) dA^ (i=x,y,z), (4) 

where t is a parameter to adjust the irregular magnetic fields strength. 

The vector potential in the point other than the lattice points are given 
by linear interpolation. The irregular magnetic fields are given by 
— ^ 

DH = rot ( A ) , (5) 

with the random vector potential A. The total fields are given as the sum 
of irregular magnetic field and average field. When we wish to simulate the 
turbulent magnetic fields with some spectrum like the Kolmogorov's one, we 
prepare a series of lattices with different lattice constants. And sum up 
the irregular magnetic fields with the weight given by the spectrum. 

The motion of the cosmic ray particles is determined by equation of 
motion; 


m 


dv 


dt 


r* 

= q(E + v x H) . 


( 6 ) 


However, in our simulation we neglect the effects from E. In other words, 
we assumed the magnetic fields is statical. This assumption can be 
justified because the velocity of cosmic ray is essentially light velocity, 
but the velocity of the magneto-hydrodynamical waves are around 10 times 
of light velocity in the interstellar circumstances. This also means that 
in our simulation, we neglect the effect of Fermi type acceleration or the 
diffusion process in energy axis. 


3 Simulation 

Our simulation is carried out as follows. We assumed that the average 
magnetic field(H^) is aline with z-axis, and the average strength of the 
magnetic field is 3 micro gauss and that of the turbulent magnetic fields 
are 1.5 micro gauss in average. (We note that what is really meaningful is 
the ratio /J)H^/Hq, which is 1/2 in our case. The absolute strength is 
needed when we translate the gyro-radius to the energy in the relation to 
the scale of the turbulence in relation to the (maximum) turbulent scale 
(=L )). We first put a number of cosmic ray particles at the origin and 
follow the motion numerically determined by the equation. The step time for 
solving (6) numerically is chosen as 1/10 of the gyro frequency and is 
changed with gyro frequency. This step time is chosen so that the error 
caused by it is small in our numerical calculation and that the simulation 
is executed within the reasonable cpu- time . Ta king th e st atistics of the 
position of cosmic ray in terms of (f± =/fy +y y and (j^ ~j(z y, we can calculate 
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K ^ from the time dependence of (jj_ and from the time dependence of (7^ 
The random magnetic field is regenerated 7 every 250 step to accelerate the 
conversion of diffusion constants to average value. 


4 Results 

In fig.l, we show a result of 
our simulation in the form of 
dimension-less quantity; K./(c x Lq) » 
where i stands for // and ]~, c right 
velocity and for maximum turbulent 
scale. The simulation of fig.l is 
carried out by assuming that the 
magnetic fields' turbulence follows 
the Kolmogorov spectrum. The crosses 
and close ones show K . . and Kj_ 
respectively. We also made a 
simulation by assuming of the only 
one turbulent scale. We found a 
little difference between them. In 
both cases, however, we can say that 
the qualitative tendencies are the 
same. Therefore the farther 
simulation is executed assuming the 
only one scale in turbulence of 
magnetic field. 



-2 -1 0 

log(R L /L 0 ) 


The whole region can be divided into several regions by the diffusion 

_i a t-i i.n n • _ • „ ; 1 4 1 r\ . „T 


tensor ^egendences on the R^. When is in the region between 10 
and 10” * xLq , the Q rgt^o 4 K^^/K_L is almost constant and bgtja are 


proportional to R 
proportional to 


tio /y/: 


.When R becomes lager than 10 


xL 0 , K 7/ is 


xL^ 


and Kj_ is L almost constant. We also gave a scale in 
fig.l to read the mfp directly frpm K,, using (1). We find in the region 

. .1 j_i ra _ t o-l 1 n ^ * ,,T ' ' f ^ 4- 1^4-4 


where the R is larger than 10 

Lf 


"'xLq',' (2) gives consistent relation 


between Kj_ and K.. assuming the mfp. However, the relation (2) is not fold 
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for R^ smaller 
where the radius is smaller th 


We ma^ 

f? 10 “ * 


see that there is an another region, 


xLq, the slope of K . , dependences on 
the R l changes to steeper (^yR^ ). However, the ambiguity or the calculation 
becomes large in this region and the reliability of simulation may be 
small . 

We also executed simulations with the condition <DH/>/Hq = 1 and 1/4 
assuming of the only one (maximum) turbulent scale. These results are shown 
in fig. 2 and in fig. 3. The simbols are saje^with fig.l^Tjae whole region is 
divided into three regions at the R ; 10 * xL« and 10 * xL^. The 


dependences of K y/ and K 


_L on R t 


v , -uw 'xLq and 10 

are the same with the case of <DH>/H 0 = 


1/2. However, there is a quantitative difference in the value of 
and in the ratio K /y /K / . They are rather rapidly varying functions of the 

<DH>. /f 


v r Kj - 


5 Summary and Discussions 

If we assume that the turbulent scale (L Q ) is 30pc and the average 
magnetic field streggjth is 3 mi^go^gauss, the energy of the cosmic ray with 
the gyro-radius 10” ’ xL is 1C; * eV for the protons. For the higher 
energy cosmic rays, the assumption of point scattering in the study of its 
propagation in Galactic space is a fairly good approximation with the 
mfp'vE . However the facter in front of the power of E is a rapidly varying 
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log(R L /L 0 ) 


log(Ri_/Lo) 


function of <DH^>. We note, however, that this is a result of the 
accumulation of small deflection by the magnetic fields and the facter 
should be obtained by different ygtjaod; for example simulation. For the 
lower energy cosmic rays than 10 eV, the diffusion constants do not 
follow the formula (l)-(3). To study the cosmic ray propagation in this 
region, we need a simulation or more sophisticated analytic methods. It is 
noted that in this region the ratio of K , , and is almost constant. If 
we assume the average strength of turbulent magnetic field is 1.5 micro 
gaijig s j the ratio is about 10. For the smaller energy cosmic ray than 
10 eV, the result of the simulation show a relatively large difference 
between different random number series. This fact might indicate that the 
other mechanism than the diffusion process like convectionshold be 
considred in the study of the propagation of cosmic ray with this energy in 
the Galaxy. 
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